Astronomy & Astrophysics manuscript no. HE2 158 paper 'final 


©ESO 2011 


September 16, 2011 





The low-metallicity QSO HE 2158-0107: A massive galaxy growing 
by the accretion of nearly pristine gas from its environment? ^ 

Bernd Husemann^, Lutz Wisotzki\ Knud Jahnke^, and Sebastian F. Sanchez^ 

^ Leibniz-Institut fur Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany 

e-mail: bhusemann@aip . de 
^ Max-Planck-Institut fur Astronomie, Konigsstuhl 17, D-69117 Heidelberg, Germany 
^ Centro Astronomico Hispano Aleman de Calar Alto (CSIC-MPIA), E-4004 Almeria, Spain 

ABSTRACT 

The metallicities of Active Galactic Nuclei are usually well above solar in their narrow-line regions, often reaching up to several 
times solar in their broad-line regions independent of redshift. Low-metallicity AGN are rare objects which have so far always been 
associated with low-mass galaxies hosting low-mass black holes (Mbh ^ lO^M©). In this paper we present integral field spectroscopy 
data of the low-redshift (z = 0.212) QSO HE 2158-0107 for which we find strong evidence for sub-solar NLR metallicities associated 
with a massive black hole (Mbh ~ 3x lO^M©). The QSO is surrounded by a large extended emission-line region reaching out to 30 kpc 
from the QSO in a tail-like geometry. We present optical and near-infrared images and investigate the properties of the host galaxy. 
The host of HE 2158-0107 is most likely a very compact bulge-dominated galaxies with a size of re ~ 1 .4 kpc. The multi-colour SED 
of the host is rather blue, indicative of a significant young age stellar population formed within the last 1 Gyr A 3cr upper limit of 
LhuigcH < 4.5 X lO^^L©^ for the H band luminosity and a corresponding stellar mass upper limit of Mbuige < 3.4 X lO^^M© show that 
the host is offset from the local black hole-bulge relations. This is independently supported by the kinematics of the gas. Although the 
stellar mass of the host galaxy is lower than expected, it cannot explain the exceptionally low metallicity of the gas. We suggest that 
the extended emission-line region and the galaxy growth are caused by the infall of nearly pristine gas from the environment of the 
QSO host. Minor mergers of low-metallicity dwarf galaxies or the theoretically predicted smooth accretion of cold (~ 10"^ K) gas are 
both potential drivers behind that process. Since the metallicity of the gas in the QSO narrow-line region is much lower than expected, 
we suspect that the external gas has already reached the galaxy centre and may even contribute to the current feeding of the black 
hole. HE 2158-0107 appears to represent a particular phase of substantial black hole and galaxy growth that can be observationally 
linked with the accretion of external material from its environment. 
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1. Introduction 

Active Galactic Nuclei (AGN) are thought to be powered by the 
accretion of ambient gas onto a supermassive Black Hole (BH) 
in the centre of their host galaxies. Because AGN can reach high 
luminosities for the most massive BHs (Mbh ~ IO^^Mq), they 
can be observed as Quasi-Stellar Objects (QSOs) up to redshifts 
beyond z>6 (Fan et al. 2003, jMortlock et al.|2011| ). The emis- 
sion lines in the Broad Line Regions (BLR) and the Narrow Line 
Regions (NLR) of AGN can be used to infer gas metallicities 
(e.g. |Hamann & Ferland| 1993 1^. QSOs are therefore considered 
as an important diagnostic for the metal enrichment and evolu- 
tion of galaxies across cosmic time (e.g. |Matteucci & Padovani| 
|T993l|Hamann & Ferland|T993l|T999l ). The metallicities in the 
BLR were found to be supersolar, reaching lOZ© in some cases, 
and to be correlated with the AGN luminosity, accretion rate, and 
BH mass (e.g. [Hamann & Fer land 1993; Shemmer et al. 2004} 
Warner et all2004t |Nagaoeral.,,2006, .Matsuoka et al.,2Qll| , 



with almost no redshift evolution. This is usually interpreted in 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max- Planck-Institut fiir Astronomie and the Instituto de Astrofisica 
de Andalucia (CSIC) and on observations collected at the European 
Organisation for Astronomical Reasearch in the Southern Hemisphere, 
Chile (program 070.B-0418). 



terms of vigorous star formation building up the massive galax- 
ies in which BHs reside and grow. 

Because the BLR is co nfined to the very circumnuclear re- 
gion on scales of <1 pc (e.g. Kaspi et al.|2QQQ ), it is not yet clear 
how the BLR metallicities are linked to the evolution of their 
host galaxies. The Narrow Line Region (NLR) , on the other 
hand, is spat ially extended on kpc scale s (e.g. 'Benn ert et al. 
2002; Schmit retaLl[2QQ3l [Bennert et a l. 2006a b; H usemann 
et al. 2008) and is certainly more representative of the host 
galaxy properties. Emission-line ratios in nearby AGN are read- 
ily observable and their values have been calibrated to metal 



abun dances v ia photoionisation models (e.g. Storchi-Bergmann 
et al. 1998; Groves et al.|2004|). Systematic studies of obscured 



AGN by [Groves et al.| |20Q6| ) and Bart h et al.[ ([2008), based 



on spectroscopic data of the Sloan Digital Sky Survey (SDSS 
[York et al.[|2000| ), revealed that AGN with NLR metallicities 
around solar and below with respect to the solar abundance set 
of [Asplund et al.[ (,2005) are extremely rare and only found in 
low-mass galaxies hosting BHs with Mbh ^ IO^Mq. Since only 
a few percent of all luminous AGN are found in ho st galax- 
ies with stellar masses less than lO^^M© (e.g. Kauff^m ann et al. 
2003 ), the BH mass-bulge mass relation Haring & Rix ( [2004 ) 
combined with the mass -metallicity relation of galaxies (e.g. 
Tremonti et al.|2Q(j4 ) thus naturally explains the rareness of low- 



metallicity AGN. 
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The QSO HE 2158-0107 at z = 0.212 was disc overed by the 
Hamburg-ESO survey (HES, Wiso tzki et alJ 20W). A radio flux 
of 1 .6 mJy at 1 .4 GHz was measured for this QSOs by the FIRST 
survey (Bec ker et al.||1995| ). In relation to its optical brightness 
(B = 16.7) the object is classified as a radio-quiet QSO. The 
SDSS targeted this QSO due to its X-ray and radio properties as 
weU and catalogued it as SDSS J220103. 13- 005300. 2. It is also 
listed in the QSO catalogue of Veron-Cetty & Veron| ( |201 0), but 
not in the catalogue of Schnei der et al.[^2010| ). The public SDSS 
spectrum of this particular QSO is of rather poor quality. In this 
paper we present evidence for HE 2158-0107 being an unusu- 
ally low-metallicity QSO for its high BH mass, with a likewise 
unusually low-mass host galaxy. We combine our own integral 
field spectrophotometry with archival optical and infrared pho- 
tometry to explore the properties of this particular QSO in more 
detail. 

Throughout this paper we assume a cosmological model with 
Ho = 70kms-^ Mpc"^ = 0.3, and = 0.7. The adopted 
cosmology corresponds to a physical scale of 3.45 kpc arc sec" ^ 
at the redshift of HE 2158-0 107 (z = 0.212). We will always 
use the AB system ( |Oke|1974| ) throughout the paper or explicitly 
refer to Vega magnitudes if required for the sake of comparison. 

2. Observations and data reduction 

2.1. Integral field spectroscopy 

We observed HE 2158-0107 with the Potsdam Multi- Aperture 
Spectro-Photometer (PMAS, |Roth et al.||2005| ) mounted on the 
3.5 m telescope of the Calar Alto Observatory in the night of 
Septe mber 6th, 2002, as part of a larger sample of radio-quiet 
QSOs ( [Husemann et al.|20'08[ Husemann et al. in prep.). Two ex- 
posures of 1800 sec each were obtained under photometric con- 
ditions at a seeing of V! 1 . 

In the lens array configuration of PMAS we covered an 8'' x 
8'' field of view centred on the QSO. 256 individual spectra were 
obtained at a spatial sampling of 0'.'5 x 0'.'5. The low resolution 
V300 grism was positioned such that the UJ3 and Ha emission 
lines could be covered in a single exposure. We determined a 
spectral resolution for this instrumental setup of 7? ~ 900 from a 
1 .15577 night sky line width of 6. 1 A FWHM. 

The reduction of the integral field data was done with the 
R3D package (Sanchez 2006]). Here we just briefly summarise the 
process as the details will be described in detail in Husemann et 
al. (in prep.) where the entire dataset will be presented. We used 
the continuum and arc lamp calibration frames taken specifically 
for the target, at the same time and airmass, to trace and ex- 
tract the 256 spectra on the CCD and to perform an absolute 
wavelength calibration. Twilight observations were obtained to 
construct a fibre flat-field to correct for the diff'erence in fibre- 
to-fibre transmission. Photometric standard star observations al- 
lowed an absolute photometric calibration of the data including 
the correction for th e atmospheric extinction at the Calar Alto 
site ( Sanchez et al.|r2007 ). A high S/N sky spectrum was ex- 
tracted from a blank sky area within the PMAS field of view, 
which we subsequently subtracted from the entire datacube. 
After correcting for diff'erential atmospheric refraction we spa- 
tially aligned the two exposures and co-added both to obtain the 
final datacube. 



2.2. Near-Infrared and optical ground-based imaging 

An H band image of HE 2158-0107 was acquired with the 
SOFI IR spectrograph and imaging camera ( [Moorwood et al. 



Table 1. Properties of the SDSS Stripe 82 images 
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1998 ) mounted on the New Technology Telescope of the La 
Silla Observatory. Ten exposures were taken in the large field 
camera mode (0^.^29 pixel" ^) consisting of 12 sub-integrations 
(NDIT) with exposure times of 5 sec (DIT) each, yielding a to- 
tal integration time of 600 sec on the source. The seeing dur- 
ing the observation was 0'.'8. We reduced these images with 
standard calibration frames using the public SOFI reduction 
pipeline provided by ESO. We measured an H band magnitude 
of 14.67 mag (Vega) for the QSO which is quite close to the 



value of 14.70 mag (Vega) reported by UKIDSS ( (Warren et al. 



2007) and indicates a good photometric calibration of the im- 



age. 



Fortunately, HE 2158-0107 is in the footprint of the 



Supern ova Survey as part of SDSS ( Ad elman-McCarthy et al. 
20081 for DR7), which covered 270 deg^ (a 2.5° wide scan along 



the celestial equator, generally called "Stripe 82") that were ob- 
served ~80 times in the autumn of the years 1998-2007. We 
stacked the multiple epoch exposures in order to reach signif- 
icantly deeper. This method was already ap plied to Strip e 82 
data for precise photometry of galaxies ( Gilba nk et al.|2010"|) and 



to study details in the structure of nearby galaxies ( Schawinski 
|et al., 2010 ). Since observations of Stripe 82 were done under var- 
ious ambient conditions, we co-added only the best 24 epochs 
with a seeing of <1V3 and a limiting sky surface brightness of 
>20 mag arcsec"^ in the r band. The combined images effec- 
tively reach ~5 times deeper (~2 mag) compared to the individ- 
ual SDSS images. Table [T] summarises the eff'ective seeing and 
limiting surface brightness for the combined images as well as 
the Galactic extinction inferred from the [Schlegel et al. (1998) 
extinction maps. 

3. Black hole mass and NLR metallicity 

3.1. QSO spectrum of HE 2158-0107 

We first obtained the nuclear spectrum of HE 2158-0107 from 
the PMAS observation using a collapsed image of the datacube 
as a spatial weighting profile for the extraction. The resulting 
QSO spectrum with a S/N of about 100 per pixel in the contin- 
uum is shown in Fig. [T] Prominent broad H/S and Ha emission 
lines are visible as well as narrow emission lines from the NLR. 
From the narrow [O iii] A5001 line ([O iii] hereafter) in the QSO 
spectrum we measured a redshift ofz = 0.2118. 

The synthetic r band magnitude of the spectrum is brighter 
by only 0.13 mag than the SDSS photometric r band magnitude, 
which may well be caused by the intrinsic variability of the QSO. 
The Fe ll complexes blue- and redward of UJ3 (e.g. Boros on &| 
Green|p? 92 ) are very weak and do not contribute significantly 
to the QSO pseudo-continuum. We find tentative evidence for 
the presence of stellar Mg and Fe absorption lines in the rest- 
frame wavelength 5 120-53 lOA (see inset in Fig. [l]). A prop- 
erly r edshifted 5 Gyr old single stellar population (SSP) spec- 
trum ( Bruzual & Chariot 2003 1 ) matches with several spectral 
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Fig. 1. Integrated QSO spectrum of HE 2158-0107 as observed with the PMAS integral field spectrograph. All prominent emission 
lines are labelled and the four telluric absorption bands in the red part of the spectrum are marked. The red and yellow curves below 
the QSO spectrum represent the arbitrarily scaled r and / band transmission curves, respectively. The inset highlights the wavelength 
region of potentially strong stellar absorption lines. The redshifted spectrum of a 5 Gyr old SSP with solar metallicity ( [Bruzual & 
|Charlot|2003 ) is overplotted as the red line in the inset for direct comparison and the main chemical species leading to the various 
absorption lines are labelled. 
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Fig. 2. Modelling of the broad RjS (left panel) 
and Ha line (right panel). The red line rep- 
resents the best-fit model to the data (black 
line) composed of individual Gaussians and 
a straight to approximate the local contin- 
uum. Individual Gaussian components for the 
B aimer lines are shown as blue lines above 
the continuum for a visual impression of 
their relative contribution. Two vertical dot- 
ted lines mark the wavelengths of the expected 
[Nil] ^^6548, 6583 lines. The residuals of the 
best-fit model are displayed in the panels be- 
low. Some residuals of the bright Na 1/15893 
night sky line are visible on top of the broad 
H/? line. Two telluric absorption bands of the 
Earth's atmosphere are clearly visible at 7600- 
7700Aand8150-8350A. 



features in this wavelength range, although some stellar absorp- 
tion lines seem to be absent. It is likely that those absorption lines 
are strongly blended in our low resolution spectrum with weak 
narrow emission lines of the QSO, like [Fevii] /15158, [Fevi] 
.15176, [Ni] .1.15177,5197, or [Fevii] .15278. Thus, the nuclear 
QSO spectrum probably contain also a small fraction of stellar 
light from the underlying host galaxy. 

3.2. Virial BH mass estimate 

BH masses of QSOs can be estimate d from their single epoch 
spectra via the so-called virial method (Peterson & Wandel 2000 ; 
[Wsterg aard 2002). The method combines the empirically de- 
rived size-luminosity rel ation for the BLR, determined via re- 
verberation m apping (e.g. Kaspi et al. 2000 Peterson et al. 2004[ 
Kaspi et al. l2005: Be ntzet al.||2006[ 



2009]), and assuming viri- 



alised motion of the BLR clouds that can be inferred from the 
broad line width. 



We used the formula also employed by |Schulze & Wisotzki 
( |20T0l ), 



Mbh = 2.57 X 10 



llOOOkms-i; \1044ergs-V ® 



(1) 



where cr^g is the line dispersion of the broad HjS line and L5100 
is the continuum luminosity at 5100 A. Equation ^ is derived 
from the empirical relation between the BLR size and the con- 
tinuum luminosity as cali brated by|Bentz et al.| ( |2009| combined 
with the prescription by |Collin et al.| ( |2()06 ) to infer the virial 
motions from cr^s adopting a virial coefficient of 3.85. By mod- 
elling the QSO spectrum with multiple Gaussians for the emis- 
sion lines and a straight line for the local continuum (Fig. [2] left 
panel), we measured an UJ3 line dispersion of cr^g = 2320 ± 
20kms"\ excluding any NLR component, and a continuum lu- 
minosity of .1L^(5100) = (8.0 ± 0.3) X 10"^^ ^-i pj-^j^ ^^ose 
measurements we estimated a BH mass of log(MBH/MQ) = 
8.6 + 0.3 for HE 2158-0 107. 

[Greene & Ho| ( |2005b] ) derived an alternative Mbh calibration 
based on the broad Hof line luminosity and width, which is less 
sensitive to host galaxy light contamination and also less suscep- 
tible to dust extinction. 



Mbh = 3.0 x 10^ 



/ FWHMHg 
\ 1000 km s"^ 



2.06 



1042 erge 



0.45 



Mo. (2) 



We modelled the Ha line with multiple Gaussian components 
(Fig. [2j right panel) and measured an Ha luminosity of (3.6 ± 
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0.4) X 10^3 erg and a FWHM of 3860 ± 70kms"^ exclud- 
ing the NLR narrow component. This yields a BH mass of 
log(MBH/Mo) = 8.4 ± 0.3 with Eq. Although the broad 
UJ3 line could be well approximated by a single Gaussian com- 
ponent, the Ha model required two Gaussian components due to 
its extended wings that are not seen for H/S. A difference in the 
line profiles is expected, because the EWHM of Uj3 is system- 
atically broader than Ha (e.g. Greene & Ho 2005b Schulze & 
|Wisotzki|2010t . 

Both BH mass estimates diff'er only by 0.2 dex and are well 
consistent with each other assuming a typical systematic uncer- 
tainty of -0.3 dex for virial Mbh estimates. This shows that a 
high-mass BH powers the QSO nucleus of HE 2158-0107. We 
adopt the average of both estimate, log(MBH/MQ) = 8.5 ± 0.3, 
for the rest of this paper. 



3.3. Evidence for an exceptionally low NLR gas-phase 
metallicity 

The estimated BH mass for HE 2158-0107 implies an expected 
bulge mass of log(Mbuige/MQ) = 11.3 ± 0.3 for its host galaxy 
adopting the empirical MsH-^buige relation of [Haring & Rix 
( [2QQ4j ). If the host galaxy has a substantial disc, the total mass 
of the system would even be higher, depending on the bulge- 
to-disc ratio. The mass-metallicity relation of galaxies implies 
that the gas-phase metallicity, usually expressed by the oxygen 
abundance, increases monotonically with stellar mass, so that 
we expect an oxygen abundance of 12 -\- log(0/H) = 9.13 ± 0.2 
from the mass-metallicity relation of [T remonti et al.| ([2004|) for 
the inferred stellar mass of HE 2158-0107. However, diff'erent 
calibrations have been used in the literature to determine the 
oxygen abundance from strong emission lines of Hii nebulae 
(e.g. 'D enicolo et "Syi[2002l [Kewley & Dopital[2QQ2l [Pettini & 



Pagel 2004n'Kobulnicky & Kewley|2004t [Tremonti et al.|2004| 



Pilyugin & Thuan|2005| ). Unfortunately, the absolute scale of the 
derived oxygen abundances depends strongly on the adopted cal- 
ibration so that the corresponding mass-metallicity relations are 
also quite diff'eren t (jLiang et al.|2006l[Kewrey & Ellison|2008[ 
Calura etal.|2009D even though the emission-line ratios are iden- 



tical. A relative comparison of abundances is therefore only pos- 
sible when they were obtained with a common calibration. 

For the radiation field of an AGN as in the case of 
HE 2158-0107, the oxygen abundance in the NLR region can- 
not be determined with the same calibrations used for Hii re- 
gions. With sta te-of-the-art photoionisation codes like CLOUDY 
( |Ferland|[T996l ), emission line ratios can be interpreted by as- 
suming a particular abundance pattern of elements, i.e. the solar 
abundance pattern, and a power-law ionising continuum of the 
AGN. Various studies reported that for the NLR the [O iii]/HyS 
ratio is primarily governed by the ionisation parameter while the 
[N ii]/Ha li ne ratio is mainly determined by the gas-phase metal- 
licity (e.g. IStorchi-Bergmann et al.]|1998t [Grov es et aLl[20Q4l 
2006| ). A complication in that interpretation is that oxygen and 



nitrogen have a diff'erent nucleosynthesis. Oxygen is predomi- 
nantly produced as a primary element in massive (young) stars 
based on the nucleosynthesis of hydrogen. Nitrogen can be pro- 
duced similarly as a primary element, but it can also have a sec- 
ondary origin in intermediate m ass (old) stars generated from the 
original metals of the star (e.g. Matteucci]|1986[ [Pilyugin et al. 
2003[[Pettini et al.|2008] ). In this simple picture the ratio of N/0 



is constant if nitrogen has a primary origin, while it depends on 
0/H if it has a secondary origin. The secondary nature of nitro- 
gen dominates at high metallic ities 12 -h log(0/H) > 8.5 (e.g. 
[Vila Costas & Edmunds[p^93|| ) where a relation between 0/H 



and N/0 can be empirically established. At the centre of mas- 
sive galaxies the metallicity of the gas is usually much higher 
than that limit. The linear relation 



log(N/0) = 0.96[12 + log(0/H)] - 9.29 



(3) 



was obtained by I Storchi-Bergmann et al. (1994) from spectro- 
scopic observations of galaxies with a nuclear starburst. 

Assuming a secondary origin of nitrogen, Storchi-Bergmann[ 



etalj(1998 ) determined a calibration for the oxygen abundances 
of the NLR by matching the [Nii]/l/l6548, 6583/HQf (= x) and 
[Oiii]/l/14960, 5007/HyS (= y) line ratios with a grid of CLOUDY 
photoionisation models described by a two-dimensional polyno- 
mial of second-order. 



12-hlog(0/H) = 8.34 -F0.212JC- 0.012^2 -0.002}; 



0.002x^3; + 6.52 X lO-y 



+221 X lO-^x/ + 8.87 X lO'^xV 
-0.11og(^e/300), 



(4) 



where the dependence on the electron density (^e) was added as 
an additional term. 

Since the oxygen abundance calibration used by 'Tremonti 
[et al.[ ( [2004| ) is also based on a match to CLOUDY photoionisa- 
tion models, but for the conditions of H II regions, the calibra- 
tion of S torchi-Bergmann et al.^ (1998) should provide oxygen 
abundances at roughly the same scale. We verify this by mea- 
suring the oxygen abundance of the NLR for all SDSS AGN 
having log 10([Oiii]/HyS) > 1.0 and log 10([Nii]/Ha^) > -0.2 
using the method of [Storchi-Bergmann et al. We find a mean 
value of 12 -h log(0/H) = 8.95 ± 0.25 which matches in abso- 
lute scale with the expected value for the mass-metallicity rela- 
tions at a stellar mass of log(M=,) ^ 10.5 where the mass di stri- 
bution of AGN host galaxy pea ks ([Kauffmann et al.[[2003| ) for 



the [Tremonti et al j (|2004| , Kobulni cky & Kewley[ ( [2004[ ), and 



Kewley & Dopita[ ([20021) calibrations whereas the other calibra 



tions are substantially off'set. We are therefore confident in using 
the mass-metallicity relation of Tremonti et al. as our reference 
for comparison with the NLR metallicity that we determine in 
the following. 

The QSO spectrum of HE 2158-0107 displays prominent 
narrow emission lines from the NLR and we included also 
Gaussian components for the narrow [Nii] /1/16548, 6583 emis- 
sion line doublet to the model of the Ha line (see Fig. [2]). 
However, no trace of the [Nii]/16583 ([Nii] hereafter) line is 
visible in the data. It is of course possible that weak [N 11] line is 
actually present and just below our detection limit considering 
the low spectral resolution and the blending with the underlying 
broad Ha line. In order to determine a realistic detection limit 
for the [N 11] line, we performed Monte-Carlo simulation as de- 
scribed in the Appendix. From those simulations we determine 
a 3cr upper limit of log([Nii]/HQf) < -1.0, or log([Nii]/HQf) < 
-1.2 at 1 cr. The location of the NLR of HE 2158-0107 in the 
log([Oiii]/HyS) vs. l og([Nii]/Ha) diagram , the so-called BPT 
diagnostic diagram ( [Baldwin et al.[p^8T] ), is shown in Fig. [3] 



The line ratios place the NLR in the domain of AGN ionisa- 
tion above the [Kewley et al.| p001[) and Kauff^mann et al.|p003| ) 
demarcation curves as expected, but at a significant offset from 
the local population of AGN towards a much lower [NiiJ/Hor 
line r atio. With the abundance calibrations of [Storchi-Bergmann 
et aL ( 1998 ) we infer an upper limit of 12 log(0/H) < 8.4 for 
the oxygen abundance in the NLR of HE 2158-0107 based on 
the [O iii]/H/5 and [N iiJ/Ha line ratios and assuming a canonical 
electron density for the NLR of ~ 1000 cm"^. This oxygen 
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Fig. 3. Standard BPT emission-line ratio diagram for 
HE 2158-0107. The measured line ratios of the NLR (filled 
red circle) and the EELR (opened red circle) are shown for the 
3cr limits on log([Nii]/HQf) and the length of the arrows are 
such that their tip corresponds to the 1 cr limit. The distribution 
of emission-line ratios for -40000 randomly selected SDSS 
galaxies are indicated by the grey scale for comparison. Triangle 
symbols indicate the NLR line ratio for a sample of AGN with 
low BH masses (Greene & Ho 2007) and blue opened squares 
represent four candidate low-metallici ty A GN in dwarf galaxies 
as reported by Izotov & Thuan (2008 ) and Izotov et al. ( 2010| ). 



abundance corresponds to Z < 0.5Zq and also agrees with the 
gas-phase metallicity of the Groves et al. (2006]) photoionisation 
models. 

Our judgement that the exceptionally low [NiiJ/Hof line ratio 
is not caused by the blending with the broad Hof line is further 
supported by line ratios of an Extended Emission Line Region 
(EELR) that we resolved around the QSO (Fig.|4]). We employed 
our software tool QDeblend^^ for the deblending of the QSO 
and extended emission using an iterative algorithm improving 
the technique introduced by [Christensen et al.| (|2006| ). Details 
on the algorithm can be found in the QDeblend"^^ manuaQ and 
its application to the PMAS IFU observation of our entire QSO 
sample will be described more detailed in Husemann et al. (in 
prep.). Briefly, a high S/N QSO spectrum from the central spaxel 
was subtracted from all other spaxels after being scaled to match 
the broad B aimer lines. Three iterations were performed to em- 
pirically decontaminate the QSO spectrum from its host con- 
tribution using a co-added spectrum of the 8 spaxels surround- 
ing the central QSO spaxel. We extracted a co-added spectrum 
from a high surface brightness region South-East of the QSO 
as indicated in Fig. |4] No broad emission lines are left in the 
extended spectrum after we applied our QSO-host deblending 
process, so that the blending of broad and narrow lines is not a 
problem for the EELR spectrum. The [N ii] doublet is still unde- 
tected in this spectrum and we determined a Scr upper limit of 
[Nii]/Ha< -0.8 ([Nii]/Ha^< -0.94 at \cr) using the same type of 
simulation method applied for the QSO spectrum. The emission- 
line ratios of the EELR are close to the unresolved NLR in the 
BPT diagram (Fig.jSj. It is likely that the AGN is the main ion- 
isation source for the extended emission as well. However, due 
to the 0.15 dex lower [Oiii]/HyS line ratio and the upper limit in 
[N ii]/HQf, we cannot strictly rule out excitation by young mas- 



sive stars from emission-line diagnostics alone. In either case, 
the limit on the [NiiJ/Hof line ratio implies that the gas-phase 
metallicity must be low over a region of several kpc, suggesting 
that the low metallicity gas is not confined to the circumnuclear 
region. 

Only very few other AGN with similar ly low oxygen abun- 
dances in the NLR gas are known to date. Groves et a l. (2006!) 
found only 40 low-metalHcity AGN among 23 000 Seyfert 2 
galaxies, all of which are in low-mass galaxies. Only two of them 
actually have gas-p hase metallici ties lower than the upper limits 
of HE 2158-0107. [Greene & H o (2007 ) constructed a sample 
of low BH mass AGN (Mbh < 2 x lO^M©) and measured the 
line ratios of the NLR which are shown in Fig. [3] for compari- 
son. They also found only 2 AGN with log(HQf/[Nii]) < -1.0 
in their NLR. Finally, 4 low-metallicit y AGN candidates with 
log( Ha/[Nii]) ~ -1.6 w ere reported by [Izotov & Thuan| ([20Q8| 
and jlzotov et al.[ ( [2010[ ), all residing in dwarf galaxies with BH 
masses in the range of the [Greene & Ho| sample. In contrast to 
the other AGN considered so far, those candidate AGN display 
emission-line ratios in the NLR that are well consistent with 
classical Hii regions powered by massive stars. Therefore, they 
should be considered as a completely diff'erent class of objects. 

HE 2158-0107 diff'ers from the above mentioned cases in 
that it is a luminous QSO with a high-mass BH. It is therefore 
the first and so far only known proper case of a low-metallicity 
QSO. Our limiting gas-phase oxygen abundance is 0.7 dex lower 
than expected for the stellar mass predicted by the BH mass- 
bulge mass relation. Either the host of HE 2158-0107 is an ex- 
treme outlier from the mass-metallicity relation of galaxies, or it 
is strongly off'set from the BH mass-bulge mass relation, or even 
both. We explore the latter aspect in the next section. 



4, The multi-colour SED of HE 2158-0107 and its 
bulge mass 

In order to directly constrain the bulge mass and to test whether 
the host galaxy of HE 2158-0107 is significantly off'set from 
the BH mass-bulge mass relation, we analysed the optical 
and infrared broad-band images to reconstruct the multi-colour 
Spectral Energy Distribution (SED) of the host galaxy. To do 
this, the emission of the QSO and the its host need to be de- 
composed, which requires a good characterisation of the Point 
Spread Function (PSF) in each band. 

Two bright stars are close to HE 2158-0107, brighter than 
the QSO by 2 mag and 1 mag in the r band, respectively. We 
empirically determined the PSF for each observation from these 
two stars. Their distances to the QSO of 6T' and 76'' already 
caused subtle PSF variations with respect to the PSF at the po- 
sition of the QSO. These variations are particularly strong for 
the SOFI image taken w ith the large field objective due to coma 
( jMoorwood et al. 1998] ). To overcome the limitations imposed 
by the field variations we analysed their PSF residuals from an- 
other star only 18'' away from the QSO, which has a similar r 
band magnitude. We used the residuals within a radius of 2" to 
make an empirical corrections to our two PSF stars. Since all the 
stars in the SDSS field are at least 0.5 mag fainter than the QSO 
in the u band, we were not able to construct a useful PSF in this 
band and excluded it from our analysis. 



4.1. Results of the QSO-host decomposition 



http://sourceforge.net/projects/qdeblend/ 



We employed GALFIT v3.0.2 ( [Peng et al.| 
to decompose the host galaxy and the point- 
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Fig. 4. [Oiii] light distribution and EELR spectrum for HE 2158-0107. Left panel: Nucleus- subtracted 40A wide [Oiii] narrow- 
band image as extracted from our PMAS datacube. The position of the QSO is highlighted by the white cross. A black box indicates 
the boundary for the extraction of a co-added spectrum. Right panel: Co-added spectrum of the EELR region as defined in the left 
panel. The expected wavelengths of various emission-lines are indicated by the vertical dotted lines. 



HE 2158-0107 by modelling the observed two-dimensional sur- 
face brightness (SB) distribution in the SDSS and SOFI broad- 
band images with an appropriate model. We considered two dif- 
ferent models. Model 1 consists of a single PSF, only to check 
whether any detectable signal from the host galaxy remains in 
the residuals even when the QSO contribution is over- subtracted. 
Model 2 is composed of a de Vaucouleur profile for the presum- 
ably bulge-dominated host and a PSF component for the QSO. 
After trying other models as well, the host turned out to be so 
faint that we fixed the SB model of the host to the de Vaucouleur 
profile to increase the robustness of the GALFIT model and the 
recovered host galaxy magnitudes (e.g. [Sanchez et al. 2004[ Kim 



|et al.||2008| ). Although other SB models for the host galaxy in- 
cluding a disc component might be possible, we emphasise that 
those models would inevitably lead to an even fainter bulge com- 
ponent. Hence, our approach provides a robust upper limit for 
the bulge luminosity. 

In Fig. [5] we present the results of the QSO-host decomposi- 
tion. In all cases, except the H band, we resolve the underlying 
host galaxy and recover the structural properties for our assumed 
host model. Interestingly, the r and / band host images reveal an 
extended tail in the light distribution, towards the South-East, 
that matches with the EELR seen in our PMAS data. We find 
that the ellipticity (e = I - b/a) and eff'ective radius of the 
best-fit de Vaucouleur profile are larger in the r and / bands than 
in the g and z bands. This is expected if line emission of the 
EELR significantly contributes to the signal in these particular 
bands (see Fig.[T]). In order to estimate the emission-line contri- 
bution of the EELR to the r and / bands, we synthesised a pure 
emission-line spectrum for the EELR based on its [O iii] flux of 
(150 ± 30) X 10"^^ erg s"^ cm"^ within an aperture of 3'' cen- 
tred on the QSO and the observed line ratios. From the synthetic 
EELR spectrum we determined an emission-line brightness in 
the r and / bands of 20.5 ± 0.2 mag and 21.6 ± 0.2 mag, respec- 
tively. This corresponds to a substantial emission-line contribu- 
tion of 18 ± 7% within the central 3'' of the r band and a rest- 
frame [Oiii] equivalent width (EW) of EW([Oiii]) = 170 + 70 A. 

Due to the apparently large nucleus-to-host ratio of > 10 in 
the z-band and the small apparent size of the host galaxy, it 
is furthermore import to estimate the systematic errors and to 
clean the structural parameters and magnitudes from systematic 
efl'ects. We performed a suite of Monte-Carlo simulations for 



various nucleus-to-host ratios and five diff'erent eff'ective radii for 
each band. We used one empirical PSF to generate 200 mock im- 
ages for a given set of parameters and analysed the image subse- 
quently with the other empirical PSF. Detailed information about 
these simulations are given in the Appendix. 

The GALFIT results of the best-fit models are listed in Table|2] 
after correcting them for Galactic extinction and for the system- 
atic biases recovered by the simulation. We determined a 3cr up- 
per limit for the host magnitude of mhost > 18.19 mag in the 
H band from our simulation. We withdraw the g band from 
this analysis because the QSO-host decomposition appears to 
be significantly aff'ected by an intrinsic PSF mismatch between 
the stars and the QSO in this band. The SEDs of the QSO and 
the stars are so diff'erent in the blue that their "eff'ective" PSFs 
are heavily aff'ected by the wavelength dependence of the see- 
ing. To quantify this eff'ect we simulated synthetic ugriz images. 
For simplicity we assumed a symmetric Gaussian profile for the 
PSF with a FWHM that is a smooth function of wavelength 
as constrained by the values in Table [T] Assuming a power- 
law spectrum for the QSO and a matched stellar spectrum for 
the stars, we constructed monochromatic images for each wave- 
length (stepsize lA). Finally, the band images were created by 
summing up the monochromatic images with its contribution set 
by the transmission curve of the band. Afterwards we scaled the 
stellar image to match with the QSO image in the brightest pixel 
and subtracted it. We find a residual flux of the order of 2% for 
the u and g bands, but only < 0.1% in the r, /, and z bands. Since 
the nucleus-to-host ratio for HE 2158-0107 is already 10 in the 
r band, such a residual flux solely due to PSF mismatch is quite 
significant. 

4.2. An SED-based stellar mass estimate 

After taking into account the systematics of the two-dimensional 
QSO-host decomposition, the emission-line contribution, and 
the foreground Galactic extinction in the diff'erent bands, 
we show the SED of the host galaxy from the rizH pho- 
tometry in Fig. |6] Even after subtracting the EELR con- 
tribution from the SED, the shape of the SED appears 
to be very blue. Scattered light of the QSO due to large 
amounts of dust within the host could significantly con- 
tribute to the apparent host light in rest-frame blue and 
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Fig. 5. Host galaxy images of HE 2158-0107 after 2D QSO-host decomposition with GALFIT. The first column represents the PSE 
(over) subtracted host images of the best-fit model including a point source only. The PSE subtracted host images and residuals of 
the best-fit model, including a point source and a de Vaucouleurs profile for the host, are shown in the second and third column, 
respectively. The radial surface brightness profiles of the data, best-fit model, and their components are presented in the fourth 
column to highlight the low contrast between nucleus and host. The orientation and position of the PMAS field of view with respect 
to the QSO is indicated by the black rectangle in the r band image and the corresponding PMAS [O iii] narrow-band image by the 
red contours. A region covering the extended tail analysed in Sect.|4.5|is marked by the blue boundary. 



UV bands (e.g. |Tadhunter et al.|1992t|Zakamska et al.|2006| ). suggests a low dust content anyway, we think that a young stellar 

We thus looked at the GALEX UV images of the Medium population as often found in bulge-dominated QSO hosts (e.g. 

Imaging Survey for HE 2158-0107. They are consistent with a [Kauffmann et al. 2003; Jah nke et al.|2004| [Sanchez et al.|2004| ) 

point source without any asymmetry that would indicate large accounts for the blue colour, 
scale scattered emission. Since the low gas-phase metallicity 
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Table 2. QSO-host decomposition results corrected for systematics 



Passband 
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[mag] 


[mag] 
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[mag] 


SDSSr 


16.18 + 0.01 


18.53 + 0.04 
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0.59 
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18.65 + 0.04 


SDSS i 


15.79 + 0.01 


18.71+0.05 


0.58 + 0.06 


0.63 


40 


18.91+0.05 


SDSSz 


15.99 + 0.02 
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0.40 + 0.10 
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24 


18.89 + 0.11 




effective wavelength [^m] 

Fig. 6. The host galaxy SED of HE 2158-0107 within a circu- 
lar aperture of y\ Open and filled symbols correspond to the 
SED with and without the contribution of the emission lines in 
the r and / band. The best-fit lOOlMyr SSP model (blue line) and 
composite 100 lMyr-i-1 1 Gyr SSP model (red line) talcen from the 
[Bmzual & Chariot, (2003 j library of spectra are shown for com- 
parison. 



In order to constrain the stellar population, w e modelled 
the SED with a grid of SSP model spectra from Bmzual & 
|Charlot| ( |2003| ) with ages in the range of 5 IVLyr to 17 Gyr and 



0.22©, 0.42©, IZ© and 2Z© metallicities. A x test was ap- 
plied to check which SSPs models were consistent with the pho- 
tometry at a 95% confidence level. We incorporated the upper 
limit in the H band by assuming a value of half the limit and a 
3cr error matching with the 3cr limit. Given the few photometric 
data points and their considerable errors we find that formally all 
SSPs with ages younger than 1 Gyr are consistent with the data. 
The best-fit SSP of solar metallicity has an age of 100 IMyr. From 
all matching SSP models we derive a minimum stellar mass of 
^buige > 3.0 X lO^lM©. This extremely low value of course just 
constrains the young stellar population and is therefore useless 
for our purposes. 

Although the light of the host galaxy is apparently dominated 
by young stars, the mass of the system could still be dominated 
by an underlying old stellar population. Thus, we determined 
how much stellar mass in an old 1 1 Gyr SSP can be added to the 
stellar mass of the young SSP, so that the composite spectrum 
still obeys the 3cr upper limit of Lbuige,i^ < 4.5 x lO^^L©^ in the 
H band. From all possible combinations we infer a maximum 
total stellar mass for the host galaxy of Mbuige < 3.4 x lO^^lM©. 



4.3. Additional mass constraints from the gas kinematics 

In Fig. [7] we present the [O iii] velocity field of the EELR which 
we derived by modelling the [O iii] doublet lines in each spaxel 





Aa [arcsec] 



offset [arcsec] 



Fig. 7. Velocity field of the [O iii] EELR around HE 2158-0107. 
Left panel: Radial velocity map inferred from the [O iii] emission 
line Doppler shift with respect to the rest- frame of the QSO. The 
black solid line represent roughly the apparent kinematic major 
axis. Right panel: Radial velocity curves extracted from the ve- 
locity map along the kinematic major axis (black open circles). 

of the nucleus- subtracted datacube with Gaussian profiles. While 
quite noisy, the overall appearance of the velocity field is consis- 
tent with ordered rotation. We constructed a synthetic longslit 
curve from the velocity map along the apparent kinematic ma- 
jor axis from all spaxels within ±0'.'5 away from it. The radial 
Doppler motion of the gas was estimated using the redshift of 
the spatially unresolved [O iii] line as the reference for the QSO 
rest- frame. From the velocity curve along the major axis we in- 
fer a maximum radial velocity of = 80 ± 20 km s"^ at a radius 
of from the nucleus. 

Assuming that the velocity field is due to a dynamically cold 
rotating gas disc within the bulge of the host galaxy, we can infer 
a rough estimate of the dynamical mass Mdyn via 



Rvi 



dyn 



G sin(/) 



(5) 



where / is the inclination of the gas disc with respect to our line- 
of-sight, R is the radius at which the radial velocity Vr is mea- 
sured, and G is the gravitational constant. The maximum of the 
radial velocity is reached somewhere between (7 kpc) and V 
(3.5 kpc), at which nearly all the mass of the bulge should be in- 
cluded considering its low eff'ective radius of re ~ 1 .4 kpc in the 
z band. Since the host galaxy has no prominent stellar disc we 
use the axis ratio of the EELR bla = 0.59 to estimate an inclina- 
tion of / = 55°. We assign a large error of ±30° to this inclination 
estimate, because the unknown geometry of the AGN ionisation 
cones with respect to the presumed gas disc can introduce a large 
systematic uncertainty to the observed hi a ratio (MulchaeyetaL] 
1996). With these very crude assumptions we estimate a possi- 
ble range in Mdyn of ~ 3 x 10^-4 x lO^^lM©, which is very well 
consistent with our SED-based upper limit on the stellar mass. 

We also tried to estimate the velocity dispersion cr^ of the 
bulge. Although we cannot directly determine cr* from the 
PIMAS spectrum of the QSO, Nelson & Whittle (1996) found 
that in Seyfert galaxies, the dispersion of the [Oiii] e mission 
line in an unresolved NLR is correlated with cr*. Thus, [Nelson] 
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Fig. 8. Left panel: MeH-^buige relation for the H band. The black triangles are measurements from Marconi & Hunt| ( [2003| ) for local 
galaxies and the dashed line is their best-fit relation. Our Scr upper limit on L//,buige for HE 2158-0107 is shown as t he red open 



their best-fit relation shown as the dashed line. An SED based 3cr upper limit for the stellar mass of the HE 2158-0107 is indicated 



circle. Middle panel: MBH-^buige mass relation. The black squares correspond to the measurements of Haring & Rix| ( [20Q4| ) with 

»7 IS 



by the red open circle. An independent estimate based on the bulge size in the z-band using the size-mass relation by |Shen et aL 
(2003 ) is shown as the blue open triangle. The light grey area represents the possible range in Mdyn of the bulge based on the ionised 
gas kinematics assuming rotational motion of a cold gas disc. Right panel: Mbh-ct* relation for local galaxies with dynamical Mbh 
measurements. Measurements by Gi iltekin et al. ( 2009 ) are shown as black round symbols and the dashed line indicate their best-fit 
relation. We used the line dispersion of the [O iii] line (green open square) in the unresolved NLR of HE 2158-0107 as a surrogate 
for (T^, following the suggestion by |Nelsoii|pOOO ). 



( |2000| ) suggested that (T[oiii] may be used as a surrogate for cr*. 
The [Oiii] line in the unresolved NLR of HE 2158-0107 can 
be described by a single Gaussian profile with a line dispersion 
of cr[oiii] = 128 ± llkms"^ after subtracting the instrumental 
resolution in quadrature. 



4.4. Is HE 2158-0107 an outlier from the BH-bulge 
relations? 

Above we collected several independent estimates for the bulge 
properties of the QSO host. Figure [5] highlights the position of 
HE 2158-0107 in comparison to the MeH-^buige,// relation of 

relation 



Marconi & Hunt' ('2003') (left panel), the MBn-M^nigc 
ot [Haring & Rix (2004), and the Mbr-o"* relation studied by 
IGtiltekin et'a l. (2009 ). We find that this QSO host is off'set from 
the local scaling relations in all three cases. 

The bulge luminosity Lbuige,// is smaller by more than 0.6 dex 
than the value indicated by the relation. One advantage of the 



Oke & Sandage, 


1968 


g. Mannucci et al. 


2001 



Chilingarian et al.|2010| ) up to z ~ 0.5 and largely independent 



of the galaxy type. We therefore do not expect that the off'set 
is caused by a systematic eff'ect due to the redshift. Since we 
adopted a robust 3(T upper limit for Lbuige,//, the true off'set from 
the relation may even be substantially larger. 

It is clear that our estimate for Mbuige from the SED is tied to 
the upper limit in Lbuige,// and thus appears similarly off'set from 
the MBH-^buige relatior]^ (middle panel of Fig. [s]). But also the 
estimated dynamical mass from the extended gas kinematics - 



^ Any stellar population model including contributions from TP- 
AGB stars ( e.g. Maraston 2005 ) will result in lower stellar masses com- 
pared to the |Bruzual & Charlot| models used here. In this case the offset 
of HE 2158-0107 from the general mass scaling relation would even be 
increased, hence our use of |Bruzual & Charlot| ( |2003| ) is conservative in 
this respect. 



while uncertain - is in agreement with the upper limit on stel- 
lar mass. Furthermore, the small apparent size of the bulge mea- 
sured in the z band supports the notion of an under-massive bulge 
with Mbuige = (1.4 ± 0.6) X IO^^Mq based on the empirical size- 
mass relation of bulge-dominated galaxies (^Shen et aL]|2003 



p007 ). It suggests that the bulge of HE 2158-0107 could be sig- 
nificantly away from MBH-Mbuige relation. Deeper high-spatial 
resolution photometry in the infrared and optical emission-line 
free bands are however required to accurately determine the SED 
of the host galaxy and to pin down the bulge luminosity and 
mass. 

The off'set in the BH mass-bulge relations appears to be most 
conspicuous in cr*. The measured value of (T[oiii] is about half of 
the expected stellar velocity dispersion {p-^ = 245 ± 44kms"^) 
from the Mbr-ct^ relation at the BH mass of HE 2158-0107 
(log(MBH/MQ) = 8.5 ± 0.3). Conversely, the measured velocity 
dispersion would correspond to a BH mass of Mbh = 2 x 10^ Mq 
adopting the Mbr-o"* relation ofl Giiltekin et al. | ([20091), smaller 
than the expected value by more than an order of magnitude. The 
discrepancy remains high even when we account for a potential 
intrinsi c off'se t Act = log((T*)-log((T[oiii]) = 0.11 dex as reported 
by Hol p009 ). However, the robustness of (T[oiii] as a surrogate 
for (T* IS strongly debated as the scatter in the relation is much 
hig her for cr ^p m] than for cr* ( Shields et al.|2003 Boroson|2005[ 
[Greene & Ho|2Q05a| ). 

Each measurement has some drawbacks and would certainly 
be considered inconclusive individually. But remarkably, all the 
independent results agree with each other and point into the same 
direction. Nevertheless, the off'set of HE 2158-0107 from the 
MBH-bulge relations is still formally consistent with the intrinsic 
scatter of the scaling relations except for the velocity dispersion. 
So far we can refer to HE 2158-0107 only as a candidate out- 
lier from the Mbh -bulge relations. The fact that the QSO spec- 
trum might be particularly suited to directly measure the stellar 
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Table 3. Photometry of the extended tail 
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Fig. 9. Gas-phase oxygen abundance against total stellar mass. 
The grey scale image correspond to a 2D histogram as con- 
structed by Tre monti et al.| p004^ from the SDSS MPA/JHU 
galaxy catalogues. The red solid line represents the correspond- 
ing best-fit mass-metallicity relation for which the dashed lines 
indicate roughly its 3cr uncertainty. Most of the AGN hosts 
among the sample of [Greene & Ho| ( [2b07| ) are within or slightly 
below this mass-metallicity relation. The blue open triangle and 
the red open circle highlight the exceptionally low metallicity of 
HE 2158-0107 with respect to the mass-metallicity relation, for 
which the stellar mass was estimated from the bulge size and the 
host SED, respectively. 



velocity dispersion from a high S/N spectrum provides a good 
opportunity to verify and pin down the off'set in the future. 

4.5. The origin of tlie low-metallicity gas 

While there is strong evidence for a significantly lower mass 
of the host galaxy than what would be expected from Mbh, its 
mass is still too high to easily explain the low metallicity of 
the gas in the NLR and EELR. The mass-metallicity relation 
decreases only by 0.2 dex from IO^^Mq to IO^^Mq in stellar 
mass as shown in Fig. [9 Thus, the gas-phase NLR metallicity 
of HE 2158-0107 is still far below the mass-metallicity rela- 
tion at the estimated stellar mass (limit) for this object. For com- 
parison we also estimated stellar ma sses and gas-phase ox ygen 
abundances for the AGN sample of Greene & Ho ( 2007 1. We 
computed rough stellar masses from the absolute g band host 
magnitudes with a mass-to-light ratio based on the g - r colour 
( Bell et al.|2003|). Since only the total g - r colour was reported 
by [Greene & Ho| ( |2007| ), but not that of the host, we assumed 
g-r ~ 0.4 as a typical value for blue- sequence galaxies ( [Blanton 
et al.|20 03]). We estimated the gas-phase oxygen abundance with 
the calibration of [Storchi-Bergmann et al. (1998 ) only when 
the line ratios were above the [Kewley et al.^ ( ,200 Ij demarca- 
tion line in the BPT diagram (cf. Fig. [3]). For classical Hii re- 
gions among the sample (i.e. objects below the'Kauff'mann et al. 
(2003 ) line) we use the abundance calibration of Pettini & Pagel 
( 2004| and converted the measured oxygen abundance to the 
sc ale of the|Tremonti et al.|C alibration following the prescription 
of |Kewley & Emson| ( |2008] ). Very few of the [Greene & Hol low 
BH mass AGN (Mbh < IO^Mq) have oxygen abundances below 
12 -h log(0/H) < 8.5 of which only three seem to have stellar 
masses comparable with that of HE 2158-0107 whereas all the 
others are consistent or only slight below the mass-metallicity 
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relation. This highlights the exceptionally low gas-phase metal- 
licity in HE 2158-0107 compared to other host galaxies of lu- 
minous AGN that have even much lower BH masses. We suspect 
that the large EELR with its intriguing tail-like geometry is the 
key to understand the unusual properties of HE 2158-0107. 

The [O III] light distribution extends beyond the PMAS field 
of view, but we can follow the structure in the Stripe 82 r 
and / band images out to 30kpc projected distance from the 
host galaxy. We determined the emission-line contribution to 
the broad-band images from a small area at the South-East part 
of the EELR (see Fig. [5]), which is still covered by the PMAS 
field of view but not contaminated by host galaxy light. We 
constructed a synthetic emission-line spectrum for that area of 
the tail including the [Oiii], HyS and Yia lines based on the ob- 
served [O III] flux, assuming an [O iii]/HyS line ratio of 7 and a 
Balmer decrement of 2.86. Furthermore, we assumed an aver- 



age [Qii]/ [Q iii] line ratio of the order of 1/3 (e.g. Villar-Martin 
et al.[[2008 ) to roughly predict the emission-line contamination 
in the g band. The surface brightness values of the tail are re- 
ported in Table [3] for both the SDSS photometry and the syn- 
thetic emission line spectrum. The corresponding SED is shown 
in Fig. [To] We find that the emission lines actually account 
for 57% of the flux in the r band whereas the Yia line ac- 
count for 25% of the / band flux. The emission lines are clearly 
dominating the r band flux, but an underlying continuum must 
also be present with a corresponding [Oiii] rest-frame EW of 
EW([Oiii]) = 1120 ± 500 A. Such a high EW for the [Oiii] line 
was only observed in the in tegrated light of 2 out of 256 type 2 
QSOs (Z akamska et al.[[2Q0 3 ). If the EELR should be powered 
by star formation, only blue compact dwarf galaxies reach such 
high equivalent widths in the [Oiii] line (e.g. Thuan & Izotov 
2005) . 

One explanation for this intriguing feature dominated by 
gaseous emission may be that it is a tidal tail due to a re- 
cent galaxy interaction. This scenario would be attractive for 
HE 2158-0107, because the oxygen abundance in interacting 
galaxies has been found to be sy s tematically lower at a given 
steflar mass ([Elfison et al.[ [20081 [Michel-Dansac et al.[[2008l 
Sol Alonso et al. 20 lO'), presumably due to inflowing low- 
metallicity gas from larger radii. Among those, the strongest di- 
lution of metallicities has been observed and theoretically pre- 
dicted for major mergers ( [Rupke et al.||2010 [Montuori et al. 
2010| ). However, while this process can eff'ectively funnel low- 
metallicity gas towards the galaxy nucleus, we find no signatures 
for a recent or ongoing major merger from the stellar light distri- 
bution at our spatial resolution of 1 kpc/pixel. Pointing to a dif- 
ferent stage of galaxy interactions, [Villar-Martin et aL ( 2010 1 re- 
cently reported the detection of a gas-rich tidal bridge of 180 kpc 
length between a type 2 QSO (z = 0.399) and a companion 
galaxy. The gas along that bridge is ionised due to ongoing star 
formation, but the QSO is capable of ionising the gas in its vicin- 
ity up to -35 kpc distance. The EW of the [O iii] emission and 
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Fig. 10. SED of the extended tail close to HE 2158-0107. The 
surface brightness of a certain region centre on the extended tail 
(cf. Fig. |5]) is shown at the effective wavelength of the corre- 
sponding passband. The initial surface brightnesses are indicated 
as the open symbols whereas the filled symbols denote the sur- 
face brightness after subtraction of the emission-line contribu- 
tion. 



the associated low gas-phase metallicity along the tidal bridge of 
that interacting system are consistent with the EELR properties 
of HE 2158-0107. However, we cannot identify a potential com- 
panion galaxy around HE 2158-0107 that could be responsible 
for such a tail. 

To put constraints on the mass of a potential interacting 
companion we estimated the ionised gas mass from the emis- 
sion line lum inosity of the entire EELR. Fo llowing the prescrip- 
tion given in [Qsterbrock & Ferland (2006]) and assuming Case 
B recombination, a low-density limit and a gas temperature of 
r 10, 000 K for ionised nebulae, the ionised gas mass is given 
by 



Mi, 



lOOcm" 



Lu/3 



1041 ergs- 



- lO^Mo 



(6) 



where Hq is the electron density of the gas. Since Hq cannot be 
estimated from our data, we assume a density of the order of 
100 cm"^ , which is a typical value for bright EELR around QSOs 
on kpc scales (e.g. Stockton et al.|2002] ). The [Oiii] luminosity 
of the entire EELR is L[Oiii] = 2.1 x lO^^ergs"^ which corre- 
sponds to an UJ3 luminosity of L^g = 3 x lO^^ergs"^ assum- 
ing an average [O iii]/Uj3 ratio of 7 throughout the entire EELR. 
Under these assumptions we estimate an ionised gas mass of 
roughly Mion ~ 3 X 10^ Mq. This is not a large amount of ionised 
gas even though cold molecular and neutral gas could increase 
the total gas mass by an order of magnitude. A tidally disrupted 
small gas-rich companion would thus be sufficient to account for 
the observed ionised gas mass. 

Considering the exceptionally low metallicity of the gas in 
the NLR of this luminous QSO, the accretion of external and 
nearly pristine gas from the environment is the most plausible 
explanation to account for the strong dilution of the gas-phase 
oxygen abundance at the galaxy centre. Our limits on the oxygen 
abundance are actually in a regime where nitrogen may be ex- 
pected to be predominantly of primary origin for which the used 
calibration of |Storchi-Bergmann et al.| is not valid any more, but 
would then even strengthen our suggestion that the inflowing gas 



is nearly pristine. Whether the infall of gas is caused by interac- 
tions with small satellite galaxies or by the "smooth" accretion 
of gas is difficult to distinguish, because our observations do not 
provide enough spatial resolution to resolve substructures in the 
EELR and are not deep enough to reliable detect any related un- 
derlying stellar continuum. 

5. Discussion and Conclusions 

Combining PMAS integral field spectroscopy and multi-colour 
broad band images, we discovered several interesting properties 
of the luminous QSO HE 2158-0107. We find strong evidence 
for a sub-solar metallicity of its NLR. Only very few AGN are 
currently known to have such low NLR metallicities, and all of 
them are low mass system with much lower BH masses than the 
value found for HE 2158-0107, log(MBH/Mo) = 8.5 + 0.3. This 
QSO clearly is a highly unusual object at least at low redshift. 

Another interesting property of this QSO is that its host 
galaxy appears to be offset from basically all known BH-bulge 
relations established for local inactive galaxies. The stellar mass 
and H band luminosity, inferred from the SED of the bulge- 
dominated host, are both at least 0.6 dex lower than expected, 
which is independently supported by the kinematics of the 
ionised gas. Clearly, the offset for a single object can always 
be explained by chance given the intrinsic scatter of the relation, 
but the low stellar mass estimated from the size-mass relation 
of galaxies and the low velocity dispersion of the nuclear gas 
suggest that HE 2158-0107 is an outlier from the local Mbh- 
bulge relations. We tentatively note that the estimated offset from 
the relation matches, by order of magnitude, with the system- 



atic offsets found at higher redshifts (e.g. [Peng et al. 


2006[|Treu 


et al.|2007 Schramm et al.|2008l|Merloni et al.|2010 


).|Jahnke & 


Maccio (20111 recently studied with numerical simulations how 



the BH mass-bulge relations may get established over cosmic 
time. They found that merging of galaxies in a ACDM universe 
alone is sufficient to produce a tight relation at z = 0, even if the 
BH and bulge masses are initially uncorrected. In fact, the local 
relation would have a much lower scatter than actually observed 
if no additional physical processes operate that increase the scat- 
ter. Thus, it is important to understand the physical mechanisms 
that lead to the observed scatter in the BH mass-bulge relations. 
It may be that HE 2158-0107 is currently offset from the rela- 
tion due to such a process or event that influences the relative 
growth of BH and bulges. 

Whether or not the low-metallicity gas in this object is 
causally connected to the higher-than-expected MBn/^buige 
tio is a key question. A striking feature of HE 2158-0107 is its 
large EELR with a tail-like geometry extending out to 30 kpc 
from the QSO. EELRs are common around radio-loud QSOs 
and are well-studied. However, it is still controversial whether 
they are produced by tidal debris from interacting galaxies, in- 
flated by galactic superwinds, or whether they may also rep- 
resent smoothly accreted cold gas from the environment (e.g. 
Stockton & MacKenty||1987[ [Fabian et al.ll 19871 [Keres et al." 
2005t [Stockton et al.||2006t |Fu & Stocktonl2009| ). In the case 
of our radio-quiet QSO, a galactic superwind can be excluded 
with high confidence, because the EELR kinematics are rather 
quiescent and a large scale radio jet is absent. A major merger 
origin of the extended tail is also unlikely considering its low 
ionised gas mass and the lack of obvious major merger signa- 
tures in the continuum. Instead we propose that the accretion of 
external low-metallicity gas explains the exceptional properties 
of this QSO and its host galaxy. This process supplies fuel for 
the growth of the bulge and presumably also of the BH. 
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The details of the presumed accretion process remain diffi- 
cult to assess with our observations. One option is that the mate- 
rial is accreted from one, possibly several low-metallicity dwarf 
galaxies in a sequence of minor merger events. If so, this pro- 
cess must have been going on sufficiently long for substantial 
amounts of low-metallicity gas to reach the NLR of the QSO. An 
alternative, but also more speculative possibility is the smooth 
accretion of gas directly from the intergalactic medium, which 
would naturally explain the observed low metallicity of the gas. 
A fraction of the accreted gas needs to be dense enough in order 
to overcome the virial shock while falling onto the galaxy and 



Birnboim 



to remain at relatively cool temperatures ~ 10^ K (e.g. 
|& Dekel 2003 ; Keres et al. 2005 ), which would then light up as 
an EELR. This process is theoretically predicted, but observa- 
tional evidence for its existence is still rare with only a few can- 
didate detections ( |Giavalisco et al . 2011; Ri baudo et al.||201l ]), 
also because the phenomenon is expected to be more common 
at high redshift than in the local universe. To shed more light on 
this issue, much deeper spectra covering the entire optical win- 
dow are needed to turn our currents limits into robust abundance 
measurements that can then be directly compared with those of 
dwarf galaxies. In this respect a study of the stellar metallicity 
of the host galaxy via absorption lines appear much more diffi- 
cult given its compactness and the blending with the QSO. On 
the other hand, the QSO provides a bright background source 
to measure abundance of different elements like O, N and C di- 
rectly from rest-frame ultraviolet metal absorption lines of the 
gas that could be compared with the abundance pattern for dwarf 
galaxies or expectations for the inter galactic medium at this low 
redshift. 

Independently of the question whether the accretion pro- 
cess is "smooth" or "clumpy", in our preferred interpretation 
HE 2158-0107 represents a particular phase of substantial 
galaxy and BH growth that can be observationally linked with 
the accretion of external gas. A fundamental difference in the 
relative BH and bulge growth might potentially occur when a 
BH is stochastically fed by an intrinsic gas reservoir within its 
host galaxy in contrast to this process where large amounts of 
nearly pristine gas are dumped onto the galaxy. More detailed 
observations are certainly needed to further constrain the prop- 
erties of this intriguing object and to elucidate on its importance 
as a "laboratory" for studying such a particular phase in the evo- 
lution of massive galaxies. 
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Fig. A.l. Fraction of 3cr [Nii] detections among 200 simulated 
mock spectra as a function of simulated HQf/[N ii] ratio. The open 
squares correspond to simulated QSO spectra and the open tri- 
angles correspond to simulated EELR spectra, respectively. The 
dashed and solid lines highlight the 95% and 68% confidence 
limits corresponding to a 3cr and Icr detection limit, respectively. 



Based on the best-fit Ha model for the QSO and EELR 
spectrum we generated mock spectra including narrow [Nii] 
/1/16548, 6583 lines, with the same profile and redshift as the 
narrow Ha line, within a certain range in Ha/[Nii]/16583 ratios 
(Ha/[N ii] hereafter) . We further added realistic Gaussian noise 
to the spectra including night sky lines using the variance vector 
matched to the observed data. 200 mock spectra were generated 
for each Hck/[N ii] ratio. 

In order to test for each mock spectrum whether the [Nii] 
lines could be recovered at a certain confidence level, we fit- 
ted the spectra with two diff'erent models, with and without the 
[Nil] components. The [Nii] lines actually add only one addi- 
tional free parameter to the model, the [Nii]/16583 flux, as the 
[Nil] .16548/[N II] .16583 flux ratio is fixed to 1/3 and the line 
width and redshift were both coupled to the narrow Ha com- 
ponent. We employed the statistical F-test to check whether the 
model including the [N ii] is a significantly better model than the 
one without based on the values for a given simulated spec- 
trum at a 95% confidence level. 



In Fig. A.l we present the fraction of spectra for which the 
model including a [N ii] component was significantly better rep- 
resenting the data as a function of the input Ha/[Nii] line ra- 
tio. At high Ha/[Nii] ratios, the [Nii] is detected in 100% of 
the spectra. This fraction drops at a certain HQf/[N ii] ratio below 
95% or 69%, which we adopt as our 3cr and cr detection limits, 
and smoothly decreasing towards zero with decreasing HQf/[N ii] 
ratio as expected. Our simulations imply a 3(T upper limit of 
Ha/[Nii] < -0.8 for the EELR spectrum (Ha/[Nii] < -0.8 Icr) 
and Ha/[Nii] < -1.0 for the QSO spectrum (Ha/[Nii] < -1.2 
Icr), respectively. 



Appendix A: Estimation of the [N ii] detection limit 
in the spectra 

The [Nil] /1/16548, 6583 doublet could neither be detected in the 
QSO nor in the EELR spectra obtained for HE 2158-0107 with 
the PMAS spectrograph. Due to the diagnostic importance of 
this line we determined robust upper limits for those lines with 
the aid of dedicated Monte Carlo simulations. 



Appendix B: Systematics of the QSO-host 
decomposition 

The two-dimensional QSO-host decomposition is a complex 
process which is not free of systematic eff'ects, especially in the 
high contrast regime. Here we present Monte Carlo simulation 
designed to investigate these systematic eff'ects specifically for 
HE 2158-0107. We infer realistic errors for each band and ob- 
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Fig.B.l. Results of Monte Carlo simulations to estimate the systematic uncertainties of the 2D decomposition with GALFIT for 
HE 2158-0107. The upper panels present the difference in simulated and observed host magnitudes (Amhost) as a function of the 
simulated host magnitude for 5 different effective radii assuming a de Vaucouleurs surface brightness profile. Detection limits (Scr) 
in host brightness for the different effective radii are indicated by the solid vertical lines in their corresponding colour. The lower 
panels show the ratio of the observed over the simulated effective radii as a function of the simulated host brightness. 



tain a robust upper limit for the host brightness in the H band. We tection limit for the host in the this band depends critically on 
created a suite of mock images for each band with a fixed QSO and corresponds to mhost,// > 18.19 for re - Qf!AA. 
magnitude corresponding to the observation. The surface bright- 
ness (SB) distribution of the underlying host galaxy was set to a 
de Vaucouleur law, the effective radii were set to 0^.^32, 0^.^44, 
0^.^55, 0^.^67 or 0'.'79, and the host magnitudes varied within a rea- 
sonable range in steps of 0.3 mag. We produced 200 mock im- 
ages for each set of parameters from the two-dimensional grid in 
parameter space. The host orientations and the axis ratios (bid) 
were randomly chosen where bja was restricted to a range of 
0.5-1.0. The simulated images were convolved with one of the 
available two PSF stars and the Gaussian noise pattern of the ob- 
served images was added before the images were analysed with 
GALFIT using the other PSF star. 

The results of the Monte-Carlo simulations, shown in 



Fig. B.l reveal that each band is subject to different system- 
atic effects. Obviously, the r band is the deepest image of the 
four with the highest S/N introducing no discernable systemat- 
ics. For the / band image, we found that increases strongly 
with the nucleus-to-host ratio, so that the was significantly 
overestimated. We suspect that this could be caused by a PSF 
mismatch between the two PSFs for this band. The recovered 
Tg in the z and H bands displays only little variation with host 
magnitude, but appears to be slightly offset from the true input 
value. 

The uncorrected host magnitudes have substantial system- 
atic errors that can only be determined with these simulations. 
Again the / band exhibits the strongest deviations from the input 
values. In all cases, except the r band, we found that Amhost is 
decreasing with increasing host brightness, which corresponds 
to a flux transfer from the PSF to the host component. The sys- 
tematic offsets and errors matched to the observed parameters of 
the host in the r, / and z bands are Amhost,r = 0.00 ± 0.02 mag, 
Amhost,/ = 0-08 ± 0.04 mag and Amhost,z = -0.17 ± 0.11, respec- 
tively. The recovered host fluxes were in all bands significantly 
above the 3cr detection limit, except in the H band. The 3(T de- 
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